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black hole quadrupole 
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In Eqs. (3.16), (3.17), (3.18), (3.24), (3.25) and (3.26) of this paper, the variable r should be replaced everywhere 
by the variable r, and the variable should be replaced everywhere by the variable 0. The definitions of f and 9 are 
given in Eq. (2.11). These replacements do not affect the any of the subsequent results in the paper. 

Also, the right hand side of Eq. (B3) is missing a term —ASL z r. 
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We analyze the effect of gravitational radiation reaction on generic orbits around a body with an 
axisymmetric mass quadrupole moment Q to linear order in Q, to the leading post-Newtonian order, 
and to linear order in the mass ratio. This system admits three constants of the motion in absence 
of radiation reaction: energy, angular momentum along the symmetry axis, and a third constant 
analogous to the Carter constant. We compute instantaneous and time-averaged rates of change 
of these three constants. For a point particle orbiting a black hole, Ryan [l5| ] has computed the 
leading order evolution of the orbit's Carter constant, which is linear in the spin. Our result, when 
combined with an interaction quadratic in the spin (the coupling of the black hole's spin to its own 
radiation reaction field), gives the next to leading order evolution. The effect of the quadrupole, like 
that of the linear spin term, is to circularize eccentric orbits and to drive the orbital plane towards 
antialignment with the symmetry axis. 

In addition we consider a system of two point masses where one body has a single mass multipole 
or current multipole of order I. To linear order in the mass ratio, to linear order in the multipole, 
and to the leading post-Newtonian order, we show that there does not exist an analog of the Carter 
constant for such a system (except for the cases of an Z = 1 current moment and an I = 2 mass 
moment). Thus, the existence of the Carter constant in Kerr depends on interaction effects between 
the different multipoles. With mild additional assumptions, this result falsifies the conjecture that 
all vacuum, axisymmetric spacetimes posess a third constant of the motion for geodesic motion. 
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I. INTRODUCTION AND SUMMARY 

The inspiral of stellar mass compact objects with 
masses /i in the range [i ~ 1 — lOOM© into massive 
black holes with masses M ~ 10 5 — 10 7 M Q is one of 
the most important sources for the future space-based 
gravitational wave detector LISA. Observing such events 
will provide a variety of information: (i) the masses and 
spins of black holes can be measured to high accuracy 
(~ 10~ 4 ); which can constrain the black hole's growth 
history [Ij; (ii) the observations will give a precise test 
of general relativity in the strong field regime and unam- 
biguously identify whether the central object is a black 
hole 0; and (hi) the measured event rate will give in- 
sight into the complex stellar dynamics in galactic nu- 
clei [l[. Analogous inspirals may also be interesting for 
the advanced stages of ground-based detectors: it has 
been estimated that advanced LIGO could detect up to 
~ 10 — 30 inspirals per year of stellar mass compact 
objects into intermediate mass black holes with masses 
M ~ 10 2 - 1O 4 M in globular clusters ||. Detect- 
ing these inspirals and extracting information from the 
datastream will require accurate models of the gravita- 
tional waveform as templates for matched filtering. For 
computing templates, we therefore need a detailed un- 
derstanding of the how radiation reaction influences the 
evolution of bound orbits around Kerr black holes 0-0] ■ 

There are three dimensionless parameters characteriz- 
ing inspirals of bodies into black holes: 

• the dimensionless spin parameter a — |S|/Af 2 of 



the black hole, where S is the spin. 

• the strength of the interaction potential e 2 = 
GM/rc 2 , i.e. the expansion parameter used in post- 
Newtonian (PN) theory. 

• the mass ratio fi/M. 

For LISA data analysis we will need waveforms that are 
accurate to all orders in a and e 2 , and to leading order 
in n/M. However, it is useful to have analytic results in 
the regimes a <C 1 and/or e 2 <C 1. Such approximate 
results can be useful as a check of numerical schemes 
that compute more accurate waveforms, for scoping out 
LISA's data analysis requirements [HHI, and for assessing 
the accuracy of the leading order in fi/M or adiabatic 
approximation [§4l0l|. There is substantial literature on 
such approximate analytic results, and in this paper we 
will extend some of these results to higher order. 

A long standing difficulty in computing the evolution 
of generic orbits has been the evolution of the orbit's 
" Carter constant" , a constant of motion which governs 
the orbital shape and inclination. A theoretical prescrip- 
tion now exists for computing Carter constant evolution 
to all orders in e and a in the adiabatic limit /i -C M 

d, EMI, but it; has not y et 

been implemented numer- 
ically In this paper we focus on computing analytically 
the evolution of the Carter constant in the regime a <C 1 , 
e <C 1, n/M <C 1, extending earlier results by Ryan 

qui. 

We next review existing analytical work on the effects 
of multipole moments on inspiral waveforms. For non- 
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spinning point masses, the phase of the I = 2 piece of 
the waveform is known to 0(e 7 ) beyond leading order 
[l6j |. while spin corrections are not known to such high 
order. To study the leading order effects of the central 
body's multipole moments on the inspiral waveform, in 
the test mass limit p -C M, one has to correct both the 
conservative and dissipative pieces of the forces on the 
bodies. For the conservative pieces, it suffices to use the 
Newtonian action for a binary with an additional multi- 
pole interaction potential. For the dissipative pieces, the 
multipole corrections to the fluxes at infinity of the con- 
served quantities can simply be added to the known PN 
point mass results. The lowest order spin-orbit coupling 
effects on the gravitational radiation were first derive d by 
Kidder [TpLthen extended by Ryan [3[l5j], Gergely [Hf, 
and Will [ijj]. Recently, the corrections of 0(e 2 ) beyond 
the leading order to the spin-orbit effects on the fluxes 
were derived [2(J Hl[ . Corrections to the waveform due 
to the quadrupole - mass monopole interaction were first 
considered by Poisson [22[ , who derived the effect on the 
time averaged energy flux for circular equatorial orbits. 
Gergely [23[ extended this work to generic orbits and 
computed the radiative instantaneous and time averaged 
rates of change of energy E, magnitude of angular mo- 
mentum |L|, and the angle k = cos (S • L) between the 
spin S and orbital angular momentum L. Instead of the 
Carter constant, Gergely identified the angular average 
of the magnitude of the orbital angular momentum, L, as 
a constant of motion. The fact that to post-2-Newtonian 
(2PN) order there is no time averaged secular evolution of 
the spin allowed Gergely to obtain expressions for L and 
k from the quadrupole formula for the evolution of the 
total angular momentum J = L + S. In a different paper, 
Gergely [l8| showed that in addition to the quadrupole, 
self-interaction spin effects also contribute at 2PN order, 
which was seen previously in the black hole perturbation 
calculations of Shibata et al. [24 |. Gergely calculated 
the effect of this interaction on the instantaneous and 
time- averaged fluxes of E and |L| but did not derive the 
evolution of the third constant of motion. 

In this paper, we will re-examine the effects of the 
quadrupole moment of the black hole and of the leading 
order spin self interaction. For a black hole, our analysis 
will thus contain all effects that are quadratic in spin to 
the leading order in e 2 and in jj,/M. Our work will extend 
earlier work by 

• Considering generic orbits. 

• Using a natural generalization of the Carter-type 
constant that can be defined for two point particles 
when one of them has a quadrupole. This facilitates 
applying our analysis to Kerr inspirals. 

• Computing instantaneous as well as time-averaged 
fluxes for all three constants of motion: energy 
E, z-component of angular momentum L z , and 
Carter-type constant K. For most purposes, only 
time-averaged fluxes are needed as only they are 



gauge invariant and physically relevant. However, 
there is one effect for which the time-averaged 
fluxes are insufficient, namely transient resonances 
that occur during an inspiral in Kerr in the vicin- 
ity of geodesies for which the radial and azimuthal 
frequencies are commensurate [Tol . |25| . The instan- 
taneous fluxes derived in this paper will be used in 
[To( | for studying the effect of these resonances on 
the gravitational wave phasing. 

We will analyze the effect of gravitational radiation re- 
action on orbits around a body with an axisymmetric 
mass quadrupole moment Q to leading order in Q, to the 
leading post-Newtonian order, and to leading order in the 
mass ratio. With these approximations the adiabatic ap- 
proximation holds: gravitational radiation reaction takes 
place over a timescale much longer than the orbital pe- 
riod, so the orbit looks geodesic on short timescales. We 
follow Ryan's method of computation 14]: First, we cal- 
culate the orbital motion in the absence of radiation re- 
action and the associated constants of motion. Next, we 
use the leading order radiation reaction accelerations that 
act on the particle (given by the Burke-Thorne formula 
[IgJ augmented by the relevant spin corrections [l4|) to 
compute the evolution of the constants of motion. In the 
adiabatic limit, the time-averaged rates of change of the 
constants of motion can be used to infer the secular or- 
bital evolution. Our results show that a mass quadrupole 
has the same qualitative effect on the evolution as spin: it 
tends to circularize eccentric orbits and drive the orbital 
plane towards antialignment with the symmetry axis of 
the quadrupole. 

The relevance of our result to point particles inspi- 
ralling into black holes is as follows. The vacuum space- 
time geometry around any stationary body is completely 
characterized by the body's mass multipole moments 
II = I ai o, ai and current multipole moments Sl — 
s ai . a 2 ...a l [27]. These moments are defined as coefficients 
in a power series expansion of the metric in the body's 
local asymptotic rest frame 28]. For nearly Newtonian 
sources, they are given by integrals over the source as 

II = I ai ,...ai = J px <ai ■ ■ -x ai> d 3 x, (1.1) 

Sl ~ ^ai,...flj — j pXpVq€pq< ai X a2 . . . x ai ^d x.(1.2) 

Here p is the mass density and v q is the velocity, and " < 
■ ■ ■ >" means "symmetrize and remove all traces". For 
axisymmetric situations, the tensor multipole moments 
Ijj (Sl) contain only a single independent component, 
conventionally denoted by /; (Si) [27|. For a Kerr black 
hole of mass M and spin S, these moments are given by 

Il + iSi = M l+1 {ia)\ (1.3) 

where a is the dimensionless spin parameter defined by 
a = |S|/M 2 . Note that Si = for even I and h = for 
odd I. 
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Consider now inspirals into an axisymmetric body 
which has some arbitrary mass and current multipoles 
Ii and Si . Then we can consider effects that are linear in 
Ii and Si for each I, effects that are quadratic in the mul- 
tipoles proportional to IiSi>, SiSf, effects that are 
cubic, etc. For a general body, all these effects can be sep- 
arated using their scalings, but for a black hole, // oc a 1 
for even I and Si oc a 1 for odd I [see Eq. (|1.3[) ]. so the ef- 
fects cannot be separated. For example, a physical effect 
that scales as 0(a 2 ) could be an effect that is quadratic 
in the spin or linear in the quadrupole; an analysis in 
Kerr cannot distinguish these two possibilities. For this 
reason, it is useful to analyze spacetimes that are more 
general than Kerr, characterized by arbitrary 7; and Si, 
as we do in this paper. For recent work on computing 
exact metrics characterized by sets of moments Z; and Si , 
see Refs. [2!| [3(| and references therein. 

The leading order effect of the black hole's multipoles 
on the inspiral is the 0(a) effect computed by Ryan [l5j |. 
This O(a) effect depends linearly on the spin Si and is 
independent of the higher multipoles Si and Ii since these 
all scale as 0(a 2 ) or smaller. In this paper we compute 
the 0(a 2 ) effect on the inspiral, which includes the lead- 
ing order linear effect of the black hole's quadrupole (lin- 
ear in I2 = Q) and the leading order spin self-interaction 
(quadratic in Si). 

We next discuss how these 0(a 2 ) effects scale with the 
post-Newtonian expansion parameter e. Consider first 
the conservative orbital dynamics. Here it is easy to see 
that fractional corrections that are linear in I2 scale as 
0(a 2 e 4 ), while those quadratic in Si scale as 0(o 2 e 6 ). 
Thus, the two types of terms cleanly separate. We com- 
pute only the leading order, 0(a 2 e 4 ), term. For the dissi- 
pative contributions to the orbital motion, however, the 
scalings are different. There are corrections to the radi- 
ation reaction acceleration whose fractional magnitudes 
are 0(a 2 e 4 ) from both types of effects linear in I2 and 
quadratic in Si. The effects quadratic in Si are due to 
the backscattering of the radiation off the piece of space- 
time curvature due to the black hole's spin. This effect 
was first pointed out by Shibata et al. (24j |. who com- 
puted the time-averaged energy flux for circular orbits 
and small inclination angles based on a PN expansion of 
black hole perturbations. Later, Gergely [l8[ analyzed 
this effect on the instantaneous and time-averaged fluxes 
of energy and magnitude of orbital angular momentum 
within the PN framework. 

The organization of this paper is as follows. In Sec. 
HU we study the conservative orbital dynamics of two 
point particles when one particle is endowed with an ax- 
isymmetric quadrupole, in the weak field regime, and to 
leading order in the mass ratio. In Sec. IIII1 we com- 
pute the radiation reaction accelerations and the instan- 
taneous and time-averaged fluxes. In order to have all 
the contributions at 0(a 2 e 4 ) for a black hole, we include 
in our computations of radiation reaction acceleration 
the interaction that is quadratic in the spin Si. The ap- 
plication to black holes in Sec. IIVI briefly discusses the 



qualitative predictions of our results and also compares 
with previous results. 

The methods used in this paper can be applied only 
to the black hole spin (as analyzed by Ryan [lj]) and 
the black hole quadrupole (as analyzed here). We show 
in Sec. [V] that for the higher order mass and current 
multipole moments taken individually, an analog of the 
Carter constant cannot be defined to the order of our 
approximations. We then show that under mild assump- 
tions, this non-existence result can be extended to exact 
spacetimes, thus falsifying the conjecture that all vac- 
uum axisymmetric spacetimes possess a third constant 
of geodesic motion. 



II. EFFECT OF AN AXISYMMETRIC MASS 
QUADRUPOLE ON THE CONSERVATIVE 
ORBITAL DYNAMICS 

Consider two point particles mi and m-i interacting in 
Newtonian gravity, where 7772 <C rni and where the mass 
mi has a quadrupole moment Qij which is axisymmetric: 



Qi 



d 3 xp(r) 



— Sin 

3 3 



< r 2 ^ 



(2.1) 
(2.2) 



For a Kerr black hole of mass M and dimensionless spin 
parameter a with spin axis along n, the quadrupole scalar 
is Q = -M 3 a 2 . 

The action describing this system, to leading order in 
m 2 /mi, is 



S 



dt 



1 2 

2^ V 



/i$(r) 



where v = r is the velocity, the potential is 
$(r) = -- - - p^Qf,, 



(2.3) 



(2.4) 



\x is the reduced mass and M the total mass of the bi- 
nary, and we are using units with G — c = 1 . We work to 
linear order in Q, to linear order in tt^/toi, and to lead- 
ing order in M/r. In this regime, the action (|2.3I) also 
describes the conservative effect of the black hole's mass 
quadrupole on bound test particles in Kerr, as discussed 
in the introduction. We shall assume that the quadrupole 
Qij is constant in time. In reality, the quadrupole will 
evolve due to torques that act to change the orientation 
of the central body. An estimate based on treating mi as 
a rigid body in the Newtonian field of 777,2 gives the scaling 
of the timescale for the quadrupole to evolve compared to 
the radiation reaction time as (see Appendix I for details) 



T C voi I mi 
m 2 



T 



M 
r 



V 

a . 



(2.5) 
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Here, we have denoted the dimensionless spin and 
quadrupole of the body by S and Q respectively, and 
the last relation applies for a Kerr black hole. Since 
n/M <C 1, the first factor in Eq. (|2.5p will be large, and 
since 1/a > 1 and for the relativistic regime M/r ~ 1, 
the evolution time is long compared to the radiation re- 
action time. Therefore we can neglect the evolution of 
the quadrupole at leading order. 

This system admits three conserved quantities, the en- 
ergy 



E = ^v 2 +/i$(r), 

the z-component of angular momentum 

L z = e z ■ (fir x v), 
and the Carter-type constant 

2Q^ 2 



(2.6) 



(2.7) 



K 



H 2 (r x v) 2 



(n-r)' 



AI 



, x2 1 2 M 



(2.8) 



(See below for a derivation of this expression for K). 



A. Conservative orbital dynamics in a 
Boyer-Lindquist-like coordinate system 

We next specialize to units where M = 1. We also 
define the rescaled conserved quantities by E = E / /.i, 
L z = L z /fi, K — K/[i 2 , and drop the tildes. These spe- 
cializations and definitions have the effect of eliminating 
all factors of \i and M from the analysis. In spherical 
polar coordinates (r, 9, ip) the constants of motion E and 
L z become 



E 



-<r 2 + r 2 8 2 + r 2 sin 2 8p 2 ) - - 
2 r 



r 2 sin 2 9<p. 



3 cos 



(2.9) 
(2.10) 



In these coordinates, the Hamilton- Jacobi equation is not 
separable, so a separation constant K cannot readily be 
derived. For this reason we switch to a different coordi- 
nate system (f , 9, ip) defined by 



Q 

r cos 9 = r cos 9 1 -\ ^ 

4r 2 



r sin 9 = f sin Oil — 



Q_ 

4r 2 



We also define a new time variable t by 



dt = 



dt. 



(2.11) 



(2.12) 



The action (|2.3[) in terms of the new variables to linear 
order in Q is 




(2.13) 



However, a difficulty is that the action (|2.13[) does not 
give the same dynamics as the original action (|2.3I) . The 
reason is that for solutions of the equations of motion for 
the action (|2.3[) . the variation of the action vanishes for 
paths with fixed endpoints for which the time interval At 
is fixed. Similarly, for solutions of the equations of motion 
for the action (|2.13l) , the variation of the action vanishes 
for paths with fixed endpoints for which the time interval 
At is fixed. The two sets of varied paths are not the 
same, since At ^ At in general. Therefore, solutions of 
the Euler-Lagrange equations for the action (|2.3p do not 
correspond to solutions of the Euler-Lagrange equations 
for the action (I2.13[) . However, in the special case of zero- 
energy motions, the extra terms in the variation of the 
action vanish. Thus, a way around this difficulty is to 
modify the original action to be 



S= dt 



1 2 
5 MV 



/it$(r) + E 



(2.14) 



This action has the same extrema as the action (|2.3[) . 
and for motion with physical energy E 7 the energy com- 
puted with this action is zero. Transforming to the new 
variables yields, to linear order in Q: 




(2.15) 



The zero-energy motions for this action coincide with the 
zero energy motions for the action (|2.14l) . We use this 
action (|2.15[) as the foundation for the remainder of our 
analysis in this section. 

The z-component of angular momentum in terms of 
the new variables (f, 8, ip, t) is 



f 2 sin 2 i 



dip 
dt 



, Q ■ 2 

1 - -zk snr < 



(2.16) 



We now transform to the Hamiltonian: 



H 



1 



9 

Pf ~ 
1 

2f2 



1 



E 



Q_ 

4f 3 



Li 



sin 



2f 4 
QE cos(2(?) 



(2.17) 
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and solve the Hamiltonian Jacobi equation. Denoting 
the separation constant by K we obtain the following 
two equations for the r and 9 motions: 



dt 



= 2E 



K 



9. 

2 



1 



2i 



and 



K 



QE cos(29). 



(2.18) 



(2.19) 



Note that the equations of motion (|2.18l) and (|2.19[) have 
the same structure as the equations of motion for Kerr 
geodesic motion. Using Eqs. ([2~T5]) . and ([2TTo) 

together with the inverse of the transformation (|2.11j) 
to linear order in Q, we obtain the expression for K in 
spherical polar coordinates: 



K 



r 4 (9 2 + sin 2 9y 2 ) + Q(r cos 9 - r9 sin 9) 2 + 



Q 



I -2 i 2h2 i 2 • 2 a • 2\ 

— — [r +r 6 +r sin Oip ) 



2Q 



cos 2 9. (2.20) 



This is equivalent to the formula (|2.8I) quoted earlier 



B. Effects linear in spin on the conservative orbital 
dynamics 

To include the linear in spin effects, we repeat Ryan's 
analysis [3, EH (he only gives the final, time averaged 
fluxes; we will also give the instantaneous fluxes). We 
can simply add these linear in spin terms to our results 
because any terms of order O(SQ) will be higher than 
the order a 2 to which we are working. The correction to 
the action (I2.3[) due to spin-orbit coupling is 



l— orbit 



= / dt 



2nSn l eijkXjik 



(2.21) 



We will restrict our analysis to the case when the unit 
vectors n, corresponding to the axisymmetric quadrupole 
Qij and to the spin Si coincide, as they do in Kerr. 

Including the spin-orbit term in the action (I2.3[) results 
in the following modified expressions for L z and K: 



and 



9 <? 

L z = n • ( M r x v) - -J [r 2 - (n • r) 2 ] , (2.22) 



K = (r x v) 2 - — n • (r x v) - ^(n • r) 2 



+Q 



( n . v )2_i v 2 

V ' 2 



(2.23) 



In terms of the Boyer-Lindquist like coordinates, the con- 
served quantities with the linear in spin terms included 
are: 



T -2-2 

L z = r sm 



§ (dfp\_2S s . n2 § _ s . n4 § 
\dt ) r 



K 



Aral 



2Q 



sin 2 6ip 2 ) -ASr sm 2 



cos 9 + Q(r cos 8 — r6 sin f 

r 



QM 



- — (r +r 



2a2 f r 2 sin 2 



The equations of motion are 



dt 



2L 2 



and 




QE cos{29). 



(2.25) 



(2.26) 



(2.27) 



III. EFFECTS LINEAR IN QUADRUPOLE AND 
QUADRATIC IN SPIN ON THE EVOLUTION OF 
THE CONSTANTS OF MOTION 

A. Evaluation of the radiation reaction force 

The relative acceleration of the two bodies can be writ- 
ten as 



a = -V$(r) + a rr , 



(3-1) 



where a rr is the radiation-reaction acceleration. Combin- 
ing this with Eqs. ([2~6]) . (pl2"2"j) and ^TM for E, L z and 
K gives the following formulae for the time derivatives of 
the conserved quantities: 



E = v a rr , 
L z = n ■ (r x a rr ), 



45 



K = 2(r x v) • (r x a rr ) n • (r x a rr ) 

r 

+2Q(n ■ v) (n ■ a rr ) - Qv ■ a rr . 



(3.2) 
(3.3) 

(3.4) 



The standard expression for the leading order radiation 
reaction acceleration acting on one of the bodies is (3~i| : 



1G 



,(6) 



32 



, 32 c(5) 

+ 4 5 £ pqb°k]p x i Vk - 



(3.5) 



Here the superscripts in parentheses indicate the number 
of time derivatives and square brackets on the indices 
denote antisymmetrization. 

The multipole moments Ijk(t) and Sjk(t) in Eq. (|3.5[) 
are the total multipole moments of the spacetime, i.e. 
approximately those of the black hole plus those due to 
the orbital motion. The expression p.5[) is formulated 
in asymptotically Cartesian mass centered (ACMC) co- 
ordinates of the system, which are displaced from the 
coordinates used in Sec. HI by an amount [28[ 



(3.6) 



This displacement contributes to the radiation reaction 
acceleration in the following ways: 
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1 . The black hole multipole moments and Si , which 
are time-independent in the coordinates used in 
Sec. HH will be displaced by <5r and thus will con- 
tribute to the (I + l)th ACMC radiative multipole 
& 

2. The constants of motion are defined in terms of the 
black hole centered coordinates used in Sec. UH so 
the acceleration a 1T we need in Eqs. (|3.2[) - Q3.4p 
is the relative acceleration. This requires calculat- 
ing the acceleration of both the black hole and the 
point mass in the ACMC coordinates using (|3. 51) . 
and then subtracting to find a rr = a^ r — a^f (14j . 
To leading order in /i, the only effect of the acceler- 
ation of the black hole is via a backreaction of the 
radiation field: the Zth black hole moments couple 
to the (I + l)th radiative moments, thus producing 
an additional contribution to the acceleration. 

For our calculations at 0(Sie 3 ), 0{I 2 e 4 ), 0(S 2 e 4 ), we 
can make the following simplifications: 

• quadrupole corrections: The fractional corrections 
linear in I 2 = Q that scale as 0(a 2 e 4 ) require only 
the effect of h on the conservative orbital dynamics 
as computed in Sec. HIK and the Burke-Thorne for- 
mula for the radiation reaction acceleration [given 
by the first term in Eq. (|3.5[) ] . 

• spin-spin corrections: As discussed in the intro- 
duction, the fractional corrections quadratic in Si 
to the conservative dynamics scale as 0(a 2 e 6 ) and 
are subleading order effects which we neglect. At 
0(a 2 e 4 ), the only effect quadratic in Si is the 
backscattering of the radiation off the spacetime 
curvature due to the spin. As discussed in item 1. 
above, the black hole's current dipole Si = Sife 
(taking the z-axis to be the symmetry axis) will 
contribute to the radiative current quadrupole an 
amount 



r>spin 3 M ry 

8 « = lI SlM 



3 3 • 



(3.7) 



quadrupole replaced by just the spin contribution 
(|3.7[) . and from Eq. (|3.8[) evaluated using only the 
orbital current quadrupole. 

With these simplifications, we replace the expression 
31) for the radiation reaction acceleration with 



16 



0(6) spin 



32 (5) S pi n 32 
45 e -' M pk XkV i A 

~f(5) orbit 



-y(5) spin 



~^ e pq{j k\ P x i Vk 



15 



S 1 S. 



i?, 



S 



13 



ni(5) spi: 
ij 



(3.9) 



To justify these approximations, consider the scaling of 
the contribution of black hole's acceleration to the orbital 
dynamics. The mass and current multipoles of the black 
hole contribute terms to the Hamiltonian that scale with 
e as 



AH ~ Sie 



21+3 



he 



21+2 



(3.10) 



Since the Newtonian energy scales as e 2 , the fractional 
correction to the orbital dynamics scale as 



AH/E ~ Sit M+l & he 



(3.11) 



To 0(e 4 ), the only radiative multipole moments that con- 
tribute to the acceleration (|3.5[) are the mass quadrupole 
-Z2, the mass octupole I3, and the current quadrupole S2 
(cf. [I3|)- Since we are focusing only on the leading or- 
der terms quadratic in spin (these can simply be added 
to the known 2PN point particle and 1.5PN linear in spin 
results) , the only terms in Eq. (I3.5[) relevant for our pur- 
poses are those given in Eq. (|3.9[) . The results from a 
computation of the fully relativistic metric perturbation 
for black hole inspirals [24[ show that quadratic in spin 
corrections to the I = 2 piece compared to the flat space 
Burke-Thorne formula first appear at 0(a 2 e 4 ), which is 
consistent with the above arguments. 



The black hole's current dipole Si will couple to 
the gravitomagnetic radiation field due to Sij as 
discussed in item 2. above, and contribute to the 
relative acceleration as 14]: 



— SlfeSy 



(5) 



(3.8) 



For our purposes of computing terms quadratic in 



the spin, we substitute S, 



spin 



for Sij 



in Eq. 



Evaluating these quadratic in spin terms requires 
only the Newtonian conservative dynamics, i.e. the 
results of Sec. ID and Eqs. (HI " (03) wit h thc 
quadrupole set to zero. 

linear in spin corrections: Contributions to these 
effects are from Eq. (I3.5[) with the current 



B. Instantaneous fluxes 

We evaluate the radiation reaction force as follows. 
The total mass and current quadrupole moment of the 
system are 



Qij — Qij ~\~ H' x i x ji 



oT ospin 

ij ij 



where from Eq. (|2.11[) 



%i€jkm'Ek£m> 



fcos8 1 



Q_ 

Af 2 



(3.12) 
(3.13) 



^ SU " V I 1 ~ 4~2 J C0S( />' fsin9 ( 1 ~ 4^2 ) SiU r- 



(3.14) 
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Only the second term in Eq. (|3.12l) contributes to the 
time derivative of the quadrupole. We differentiate five 
times by using 



d 



1 + ^008(20) 



d_ 

dV 



(3.15) 



to the order we are working as discussed above. Af- 
ter each differentiation, we eliminate any occurrences of 
dip/dt using Eq. (|2 . 24[) . and we eliminate any occurrences 
of the second order time derivatives d 2 f/dP and d 2 8/dP 



in favor of first order time derivatives using (the time 
derivatives of) Eqs. (|2.26|) and (|2.27|) . For computing the 
terms linear and quadratic in Si, we set the quadrupole 
Q to zero in all the formulae. We insert the resulting ex- 
pression into the formula (13.91) for the self-acceleration, 
and then into Eqs. ([3~2]) - (|3~4|) . We eliminate (df/di) 2 , 
(d9/dt) 2 , and (dip/dt) in favor of A, L z , and A using 
Eqs. (|2.24[) - (|2.27|) . In the final expressions for the in- 
stantaneous fluxes, we keep only terms that are of O(S), 
0(Q) and 0(S 2 ) and obtain the following results: 
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and 
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C. Alternative set of constants of the motion 

A body in a generic bound orbit in Kerr traces an 
open ellipse precessing about the hole's spin axis. For 
stable orbits the motion is confined to a toroidal region 
whose shape is determined by E, L z , K. The motion 
can equivalently be characterized by the set of constants 
inclination angle I, eccentricity e, and semi-latus rectum 
p defined by Hughes (32j. The constants L, p and e are 
defined by cos i = L z /vK, and by r± = p/(l ± e), where 
f± are the turning points of the radial motion, and r 
is the Boyer-Lindquist radial coordinate. This param- 
eterization has a simple physical interpretation: in the 
Newtonian limit of large p, the orbit of the particle is an 
ellipse of eccentricity e and semilatus rectum p on a plane 
whose inclination angle to the hole's equatorial plane is 
i. In the relativistic regime p ~ M, this interpretation of 
the constants e, p, and i is no longer valid because the 
orbit is not an ellipse and i is not the angle at which the 
object crosses the equatorial plane (see Ryan 14] for a 
discussion). 

We adopt here analogous definitions of constants of 
motion l, e and p, namely 



cos(i) = L Z /VK, 
P 



lie 



r±. 



(3.19) 
(3.20) 



Here K is the conserved quantity (|2 .23[) or (|2 . 25[) , and f± 
are the turning points of the radial motion using the f 



coordinate defined by Eq. (|2.11l) . given by the vanishing 
of the right-hand side of Eq. (j2~26l) . 

We now rewrite our results in terms of the new con- 
stants of the motion e, p and l. We can use Eq. (|2.26[) 
together with the equations (|3.19[) and (I3.20[) to write E, 
L z and K as functions of p, e and i. To leading order in 
Q and <S* we obtain 
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(3.21) 

, (3.22) 
Q cos 2 i 

2 

P 

(3.23) 



As discussed in the introduction, the effects quadratic in 
S on the conservative dynamics scale as 0(a 2 e 6 ) and thus 
are not included in this analysis to 0(a 2 e 4 ). 

Inserting these relations into the expressions (|3. 161) 
(|3~T8|) gives, dropping terms of O(QS), 0(Q 2 ) and 
0(QS 2 ): 
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D. Time averaged fluxes 

In this section we will compute the infinite time- 
averages (E), (L z ) and (K) of the fluxes. These averages 
are defined by 



(E) 



1 f T/2 ■ 
lim - / E{t)dt. 



T/2 



(3.27) 



These time-averaged fluxes are sufficient to evolve or- 
bits in the adiabatic regime (except for the effect of res- 



onances) [12], |25j. In Appendix II, we present two dif- 
ferent ways of computing the time averages. The first 
approach is based on decoupling the f and motion us- 
ing the analog of the Mino time parameter for geodesic 
motion in Kerr [l2j |. The second approach uses the ex- 
plicit Newtonian parameterization of the orbital motion. 
Both averaging methods give the following results: 
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Using Eqs. (|3 . 2 1 1) and (|3.23[) . we obtain from (|3.28l) - (I3.30|) the following time averaged rates of change of the 
orbital elements e, p, t: 
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IV. APPLICATION TO BLACK HOLES 



B. Comparison with previous results 



A. Qualitative discussion of results 



The above results for the fluxes, Eqs. (|33T|) . (|3~32l 
and (|3.33p show that the correction terms at 0(a 2 e 4 ) 
due to the quadrupole have the same type of effect on the 
evolution as the linear spin correction computed by Ryan: 
they tend to circularize eccentric orbits and change the 
angle i such as to become antialigned with the symmetry 
axis of the quadrupole. 

The effects of the terms quadratic in spin are quali- 
tatively different. In the expression p. 281) for (E), the 
coefficient of cos(2t) due to the spin self-interaction has 
the opposite sign to the quadrupole term, while the terms 
not involving i have the same sign. The terms involving 
cos(2t) in Eq. flOO) for (K) of 0(Q) and 0{S 2 ) terms 
have the same sign, while the terms not involving i have 
the opposite sign. The fractional spin-spin correction to 
(L g ), Eq. (|3.29j) . has no i-dependence, and in expres- 
sion p. 331) for (i), the dependence on l of the two effects 
0(Q) and 0(S 2 ) is different, too. This is not surprising 
as the 0(Q) effects included here are corrections to the 
conservative orbital dynamics, while the effects of 0(S 2 ) 
that we included are due to radiation reaction. 



The terms linear in the spin in our results for the time 
averaged fluxes, Eqs. (|3.28j) - (|3.33j) . agree with those 
computed by Ryan, Eqs. (14a) - (15c) of (IE} , and with 
those given in Eqs. (2.5) - (2.7) of Ref. (33[, when we use 
the transformations to the variables used by Ryan given 
in Eqs. (2.3) - (2.4) in 0). 

Equation (|3.28j) for the time averaged energy flux 
agrees with Eq. (3.10) of Gergely [H| and Eq. '(4.15) 
of [18| when we use the following transformations: 
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where A, L, k, 6, ipo and ipi are the quantities used by 
Gergely. The first relation here is obtained from the turn- 
ing points of the radial motion as follows. We compute 
r± in terms of E and K and map these expressions back 
to r using Eqs. (|2.11[) . The result can then be com- 
pared with the turning points in Gergely's variables, Eq. 
(2.19) of [23|], using the fact that E is the same in both 
cases. Instead of the evolution of the constants of motion 
K and L z , Gergely computes the rates of change of the 
magnitude L of the orbital angular momentum and of the 
angle k defined by cosk = (L • S)/L. Using the trans- 
formations (|4.1[) - (|4.4p and the definition of n we verify 
that our Eq. (|3.29[) agrees with the ( L 7 ) computed using 
Gergely's Eqs. (3.23) and (3.35) in Q and Eq. (4.30) 
of O. 

In the limit of the circular equatorial orbits analyzed 
by Poisson [22j |. our Eq. (I3.28|) agrees with Poisson's Eq. 
(22) when we use the transformations and specializations: 
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(4.6) 
(4.7) 
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where v and a a are the variables used by Poisson and the 
relation (|4.5I) is obtained by comparing the expressions 
for the constants of motion in the two sets of variables. 

The main improvement of our analysis over Gergely's 
is that we express the results in terms of the Carter-type 
constant K, which facilitates comparing our results with 
other analyses of black hole inspirals. Our computations 
also include the spin curvature scattering effects for all 
three constants of motion; Gergely [l8[ only considers 
these effects for two of them: the energy and magnitude 
of angular momentum, not for the third conserved quan- 
tity. 

When we expand Eq. (13.28)) for small inclination an- 
gles and specialize to circular orbits, then after converting 
p to the parameter v using Eq. (|4.5[) . we obtain 



32 
5p 5 
32 
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2p 2 



48 
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This result agrees with the terms at 0(a 2 v 4 ) of Eq. (3.13) 
of Shibata et al. (24|, whose calculations were based on 
the fully relativistic expressions. This agreement is a 
check that we have taken into account all the contribu- 
tions at 0(a 2 e 4 ). The analysis in Ref. [24| could not dis- 
tinguish between effects due to the quadrupole and those 
due curvature scattering, but we can see from Eq. (|4.9j) 
that those two interactions have the opposite dependence 
on 1. Comparing (|4.9[) with Eq. (3.7) of [24{ (which gives 
the fluxes into the different modes (I = 2, m, n), where m 
and n are the multiples of the (p and 9 frequencies), we see 
that the terms in the (2, ±2, 0) and the (2, ±1, ±1) modes 



are entirely due to the quadrupole, while the spin-spin in- 
teraction effects are fully contained in the (2, ±1, 0) and 
(2,0,±1) modes. 



V. NON-EXISTENCE OF A CARTER-TYPE 
CONSTANT FOR HIGHER MULTIPOLES 

In this section, we show that for a single axisymmctric 
multipole interaction, it is not possible to find an ana- 
log of the Carter constant (a conserved quantity which 
does not correspond to a symmetry of the Lagrangian), 
except for the cases of spin (treated by Ryan [15J) and 
mass quadrupole moment (treated in this paper). Our 
proof is valid only in the approximations in which we 
work - expanding to linear order in the mass ratio, to 
the leading post-Newtonian order, and to linear order in 
the multipole. However we will show below that with 
very mild additional smoothness assumptions, our non- 
existence result extends to exact geodesic motion in exact 
vacuum spacetimes. 

We start in Sec. IV Al by showing that there is no co- 
ordinate system in which the Hamilton- Jacobi equation 
is separable. Now separability of the Hamilton- Jacobi 
equation is a sufficient but not a necessary condition for 
the existence of a additional conserved quantity. Hence, 
this result does not yield information about the existence 
or non-existence of an additional constant. Nevertheless 
we find it to be a suggestive result. Our actual derivation 
of the non-existence is based on Poisson bracket compu- 
tations, and is given in Sec. IV Bl 



A. Separability analysis 

Consider a binary of two point masses mi and m2, 
where the mass mi is endowed with a single axisymmet- 
ric current multipole moment Si or axisymmetric mass 
multipole moment In this section, we show that the 
Hamilton- Jacobi equation for this motion, to linear order 
in the multipoles, to linear order in the mass ratio and to 
the leading post-Newtonian order, is separable only for 
the cases Si and I2. 

We choose the symmetry axis to be the z-axis and write 
the action for a general multipole as 



S= dt 



+ r 2 sin 2 9ip 2 I + - 



■f(r,e) + g(r,6)<p + E] 



(5.1) 



For mass moments, g(r,6) = 0, while for current mo- 
ments f(r,9) — 0. For an axisymmetric multipole of 
order Z, the functions / and g will be of the form 



, (r ,„_2Mg-fi, 9( „ 



diSi sin 9dgPi (cos 9) 



(5.2) 

where P; (cos 9) are the Legendre polynomials and C\ and 
di are constants. We will work to linear order in / and g. 
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In Eq. (|5.1[) . we have added the energy term needed when 
doing a change of time variables, cf. the discussion before 
Eq. (|2.14|) in section ITTTl Since ip is a cyclic coordinate, 
p v = L z is a constant of motion and the system has 
effectively only two degrees of freedom. Note that in the 
case of a current moment, there will be correction term 
in L,: 



L z = r 2 sin 2 9ip + g (r,6). 



(5.3) 



Next, we switch to a different coordinate system 
(f, 9, ip) defined by 



r = f + a(r,9,L z ), 
9 = 9 + /3{r,9,L z ), 



(5.4) 
(5.5) 



where the functions a and /3 are yet undetermined. We 
also define a new time variable t by 



dt 



l + 7(f,0,L z ) dt 



(5.6) 



where we have denoted 



Ci 

c 2 
c 3 

V 



J(r, 9) [1 + 7 - 2a, r ] = 1 + 7 - 2a, r + j, (5.9) 
J(r,8) 



l-^-2/3. e + 7 
r 



= 1 - 



2a 



- 2/3, e + 7 + j, 



(5.10) 



J(r, 9) [-a 6 - r 2 /3 r ] = -a fi - r 2 /3 r , (5.11) 
Li 2a 



J(r,6) 



2r 2 sin 2 



;(l + 7 



2/3 cot 9) 



1 

r 

f 



(l--+ 7 )-£(l + 7) 



gLz 



r z sin 



L- 



2r 2 sin 



:(l + 7 



2a 
r 



2/3 cot 9 + j) 



-£(l+7+j)- -(1- -+7 + i) 
r r 



r 2 sin 2 1 



(5.12) 



The unperturbed problem is separable, so make the 
perturbed problem separable, we have multiplied the 
Hamilton- Jacobi equation by an arbitrary function 
J(r,9), which can be expanded as J(r,9) = 1 +j(r,9), 
where j (r, 9) is a small perturbation. 

To find a solution of the form W = W r (r) + We (9) , we 
first specialize to the case where C3 = 0: 



Since we work to linear order in / and g, we can work 
to linear order in a, /3, and 7. We then compute the 
action in the new coordinates and drop the tildes. The 
Hamiltonian is given by 



H 



lp 2 (l + 7-2a r ) + ||(l 



P, 



1 (-«, 



r 

T 2 

2r 2 sin 2 9 



(1 + 7 



-V-+7)-/ 
r r 



2a 
r 

r 2 /3 r ) - £7(1 + 7) 

2a 
r 

gL z 



2/3 e+ 7) 



2,3 cot 6>) 



r 2 sin 2 t 



(5.7) 



and the corresponding Hamilton- Jacobi equation is 



= 



+2 



fdW\ (dW_\ C3 



2V, 



(5.8) 



- C 3 = P, r 



a a = 0. 



(5.13) 



We differentiate Eq. f|5.8[) with respect to 0, using Eq. 
(f5TH|) to write (dWr/dr) 2 in terms of (dWe/<i(9) 2 and then 
differentiate the result with respect to r to obtain 



fdWeY 
V d9 j 

+2d r 



Or 



dgC2 dgCi 



C2 C\ 

d e V r 2 VdgC x 



Co 



C\C'2 



(5.14) 



Expanding Eq. (|5 . 14[) to linear order in the small quan- 
tities then yields the two conditions for the kinetic and 
the potential part of the Hamiltonian to be separable: 



d r de [ 2a 
L 2 



2a 
r 



sin 



(2/3 r cot 2 



2/3 e J , (5.15) 
- 3/3^0 cot 9 + P <r esc 2 9) 



Lz 



sin ( 
-d r d g 



-d r 



a,e 
r 



a. r e 



-d r 



r ( 2a 



T Pz(cos( 



diSiL z 



r' sm ( 



de Pi (cost 



2Er 2 a r6 



(5.16) 



15 



where we have used Eq. (I5.2|) for / and g. Therefore, the 
following conditions must be satisfied: 

a 

MAO) — N(r) = -+/3 e ~2a r , (5.17) 
r 

Mi(0) = 2f3 cot 2 + (3 esc 2 + /3 ee 

-3/3, e cot 0, (5.18) 

M 2 (0) = r 2 <9 r (r 2 /3 r ), (5.19) 

M 3 (0) = 2ra,re-a,e + -^r%P/(cos0) 



<9 (csc0<9 e P(cos0)). (5.20) 



Here, the functions M and A~ are arbitrary integration 
constants. 

Solving the condition for the kinetic term to be sep- 
arable, Eq. (|5.1?p . together with Eq. (|5. 13[) gives the 
general solution that goes to zero at large r as 



.4 



Y cos(n0 + v), 



A 



(3 = — -sin(n0 + z/), 



(5.21) 
(5.22) 



where A and v are arbitrary and n is an integer. These 
functions must satisfy the conditions (|5 . 1 8[) - (|5.20[) in 
order for the potential term to be separable as well. To 
see when this will be the case, we start by considering Eq. 
(|5.20j) . Substituting the general ansatz a = ai(r)a,2(0) 
shows that a' 2 — Pi or a' 2 = (csc0 P[)' depending on 
whether a mass or a current multipole is present. The 
function a\(r) is then determined from 



= 2ra' 1 — ai 



(l-i) 



| cili/i 
ydiSiL z /r 



Hence, 



a i 



[c l I l /{2l)]r^ 
[diSiL z /(2l + l)]r- 1 



so that we obtain for mass moments 



a 



cjJi Pi (cost 



dhPl (cos < 



21 r 1 - 1 ' " 2l 2 
and for current moments 

diSiL z csc9P{(cos9) 



P = 



21 + 1 r l 

diSiL z (csc0P/(cos0))' 



(2I + 1)(Z + 1) 



r i+i 



(5.23) 

(5.24) 

(5.25) 

(5.26) 
(5.27) 



where we have used the condition (|5 . 1 3[) to solve for /?. 

Substituting this in Eq. (j5.19[) determines that I = 2 
for mass moments and 1 + 1 = 2 for current moments. For 
an I = 2 mass moment, conditions (I5.17|) and f|5 . 1 8[) are 
satisfied as well, with n = 2 and v = 0. For the case of an 
I = 1 current moment, the extra term in H is independent 



of anyway. But for any other multipole interaction, 
the Hamilton- Jacobi equation will not be separable. For 
example, for the current octupole Sijk, the last term in 
Eq. (I5.7[) is proportional to 5 l 3 L z (5cos 2 - l)/r 5 and 
is therefore not separable. From Eq. (|5.2I) one can see 
that, for a general multipole, the functions / or g contain 
different powers of cos appearing with the same power 
of r since the Legendre polynomials can be expanded as 



P;(C0S( 



N 

T> o 2i n \(l-n)\(l 



(-l) n (2l-2n)\ 
2nV. 



(cos0) 



l-2n 



(5.28) 



where N = 1/2 for even I and N = (I + l)/2 for odd I. 
It will not be possible to cancel all of these terms with 
(15~2H - (I5~22|) for I > 2. 

The case when C3 is non- vanishing will only be sepa- 
rable if all the coefficients are functions of r or of only, 
and if in addition, the potential also depends only on r or 
on 0. Achieving this for our problem will not be possible 
because the potential cannot be transformed to the form 
required for separability. 



B. Derivation of non-existence of additional 
constants of the motion 

In this subsection, we show using Poisson brackets that 
for a single axisymmetric multipole interaction, to linear 
order in the multipole and the mass ratio, a first integral 
analogous to the Carter constant does not exist, except 
for the cases of mass quadrupole and spin. 

Suppose that such a constant does exist. We write the 
Hamiltonian corresponding to the action (|5.1[) as H = 
Hq + SH and the Carter-type constant as K = Kq + 
5K(p r ,pe,L z ,r,6), where 



H 
SH 
K 



9 
Pe 



2r 2 



2r 2 sin 2 

all „ , „. diSiL 
_^. Pi(cos ^_ ____ 



r l+2 gm Q 



sin 2 



(5.29) 
d e Pi (cos 0) (5.30) 
(5.31) 



Computing the Poisson bracket gives, to linear order in 
the perturbations 

= {H ,5K} + {SH,K } (5.32a) 
= j t SK + {SH,K }, (5.32b) 

where we have used that {Hq, Kq} = and the fact that 
{Hq,SK} = d(5K)/dt. Here, d/dt denotes the total time 
derivative along an orbit (r(t),9(t),p r {t),pg(t)) of H in 
phase space. The partial differential equation (|5.32al) for 
SK thus reduces to a set of ordinary differential equa- 
tions that can be integrated along the individual orbits 
in phase space. 
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The unperturbed motion for a bound orbit is in a 
plane, so we can switch from spherical to plane polar co- 
ordinates (r,ip). In terms of these coordinates, we have 
H a = p?/2+p 2 b /2, K = and cos 6 = sintsin(^+V>o), 
with cost = L z j\f~K and the constant ipo denoting the 
angle between the direction of the periastron and the 
intersection between the orbital and equatorial plane. 
Then Eq. becomes 



±6K 

dt 



2p^ diSiL z 



ri + 2 



(5.33) 

c\,P(sinisin(-0(t) + ip )) 
cos(^(t) + ipo) 



-^^^(sintsinWt) + ^d)). (5.34) 



For unbound orbits, one can always integrate Eq. 
(I5.33[) to determine SK. However, for bound periodic 
orbits there is a possible obstruction: the solution for 
the conserved quantity Kq + SK will be single valued if 



and only if the integral of the source over the closed orbit 
vanishes, 



Torb 



■q(t)dt = 0. 



(5.35) 



Here, T or b is the orbital period. In other words, the par- 
tial differential equation (|5.32l) has a solution SK if and 
only if the condition (|5.35l) is satisfied. This is the same 
condition as obtained by the Poincare-Mel'nikov-Arnold 
method, a technique for showing the non-integrability 
and existence of chaos in certain classes of perturbed dy- 
namical systems [35j . 

Thus, it suffices to show that the condition (|5.35|) is 
violated for all multipoles other than the spin and mass 
quadrupole. To perform the integral in Eq. (|5.35l) . we use 
the parameterization for the unperturbed motion, r — 
K/(l + ecosip) and dt/dip = K 3 / 2 /(l + e cos-0) 2 , so that 
the condition for the existence of a conserved quantity 
K + SK becomes 



/.2tt 



dip 



T . ,\l-ia vi- ■ / / i i \\ diSiL z , fd^Pi {sin Lsm(ifi + ip )) s 

ciIAl + ecostp) o^PiismLSimtp + ip )j 7 (1 + ecostp) 1 

K sin l 



cos(tp + ip ) 



In terms of the variable x = V 1 + ~ 7r/2, Eq. (|5.36p can be written as 

f 2jT d 
= / d-xcih [1 + e(sin'0o cosx — cost/'o sin%)l — — P; (sin £ cosx) 
Jo ax 

f 2 \ diSiLz . , d ( 1 d 

+ / "X — : 1 + e(sm«o cos y — cos Wo sin y) -r- — — P;(sintcos y) 

Jo sin 4 dx \smxdx 



= 0. 
(5.36) 



(5.37) 



Inserting the expansion (|5.28l) for p(cosx), taking the derivatives, and using the binomial expansion for the first term 
in Eq. (jOT) . we get 



N l-l 



= Q/ i ^^A i „ jfc e J (sint)^ 2 "(sinV'o) fc (cosV'o)^ fe / d X (si 



diSiL 
K 



n=0 j=0 

N I 



sin x 



Y^B lnjk e^'(sinO'- 2n - 1 (sin^ ) fc (cos^o)^ fe f " d X (sin X )^ fe+1 (cos x)^'- 3 "" 2 



ij'-k+l. 



cosx) 



fc+Z-2n-l 



(5.38) 



n=0 j=0 

The coefficients A;„fcj and are 

1)1(2/ -2n)! 



A-lnkj — 



2 l n\{l - 1 - j)lk\(j - fc)!(Z - n)!(Z - 2n - 1)! : 



5, 



(rife J 



(_1)"+*7!(2Z - 2n)! 



2 z n!(Z - j)!fc!(j - fe)!(Z - rc)!(Z - 2n - 2)! ' 



(5.39) 



The only non- vanishing contribution to the integrals in Eq. f|5 . 38[) will come from terms with even powers of both 
cos x and sin x- These can be evaluated as multiples of the beta function: 



N l-l 



= C/7/ ^2^2C lnjk e J (smi) 1 2n (sin-0 o )' £ (cOS7/>o)' J * %-fc+l),cvcn %+fe-l),eve D 
n=0 j=0 

N I 



diSiL z 
K 



^2^2 Di njk e j (sin l) 1 2,1 ^sin Vo^cos V>o) J fc %-fe+i) 



e\ < ii + ,evei i 



(5.40) 



n=0 j=0 
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Here, the coefficients are 

_ 2r(i-| + i)r(f + i-n) _ 2r(§-| + i)r(f + 

c ^ fc - r(| + i-n + i) ^ ln3k ~ r(i + i-n + |) ( J 



Eq. (|5.40p shows that for even Z, terms with j =even 
(odd) and k =odd (even) give a non-vanishing contribu- 
tion for the case of a mass (current) multipole, and hence 
Kq+SK is not a conserved quantity for the perturbed mo- 
tion. Note that terms with j =even and k =odd for even 
I occur only for I > 3, so for I — 2 the mass quadrupole 
term in Eq. (|5.40j) vanishes and therefore there exists an 
analog of the Carter constant, which is consistent with 
our results of Sec. |TT] and our separability analysis. For 
odd I, terms with j =odd (even) and k =even (odd) are 
finite for li (Si). Note that for the case I = 1 of the spin, 
the derivatives with respect to \ in Eq. (|5.37|) evaluate to 
zero, so in this case there also exists a Carter-type con- 
stant. These results show that for a general multipole 
other than I2 and Si, there will not be a Carter- type 
constant for such a system. 



1. Exact vacuum spacetimes 

Our result on the non-existence of a Carter-type con- 
stant can be extended, with mild smoothness assump- 
tions, to falsify the conjecture that all exact, axisymmet- 
ric vacuum spacetimes posess a third constant of the mo- 
tion for geodesic motion. Specifically, we fix a multipole 
order I, and we assume: 

• There exists a one parameter family 

(M,g ab (X)) 

of spacetimes, which is smooth in the parameter A, 
such that A = is Schwarzschild, and each space- 
time gab (A) is stationary and axisymmetric with 
commuting Killing fields d/dt and d/d<fi, and such 
that all the mass and current multipole moments of 
the spacetime vanish except for the one of order I. 
On physical grounds, one expects a one parameter 
family of metrics with these properties to exist. 

• We denote by -ff(A) the Hamiltonian on the tan- 
gent bundle over M for geodesic motion in the met- 
ric gab(A). By hypothesis, there exists for each A 
a conserved quantity M(A) which is functionally 
independent of the conserved energy and angular 
momentum. Our second assumption is that M (A) 
is differentiable in A at A = 0. One would expect 
this to be true on physical grounds. 

• We assume that the conserved quantity M(A) is 
invariant under the symmetries of the system: 



where £ and if are the natural extensions to the 8 
dimensional phase space of the Killing vectors d/dt 
and d/d(j). This is a very natural assumption. 

These assumptions, when combined with our result of 
the previous section, lead to a contradiction, showing 
that the conjecture is false under our assumptions. 

To prove this, we start by noting that M(0) is a con- 
served quantity for geodesic motion in Schwarzschild, so 
it must be possible to express it as some function / of 
the three independent conserved quantities: 



M(0) = f(E,L z ,K o 



(5.42) 



Here E is the energy, L z is the angular momentum, and 
K a is the Carter constant. Differentiating the exact re- 
lation {H(X), A/(A)} = and evaluating at A = gives 

{H ,Mi} = ^ { E,Hi} + §- { L z ,Hi} + -^- { K 0> Hi}, 

(5.43) 

where H Q = H(0), H 1 = H'(0), and Mi = M'(0). As 
before, we can regard this is a partial differential equa- 
tion that determines Mi, and a necessary condition for 
solutions to exist and be single valued is that the integral 
of the right hand side over any closed orbit must vanish: 



^{E, Hi} + §^{L z ,Hi} + ^{Ko, Hi} 



0. 



dK 

(5.44) 

Now strictly speaking, there are no closed orbits in 
the eight dimensional phase space. However, the ar- 
gument of the previous section applies to orbits which 
are closed in the four dimensional space with coordinates 
(r,9,p r ,pg), since by the third assumption above every- 
thing is independent of t and </>, and pt and p$ are con- 
served. Here (t, r, 9, <fi) are Schwarzschild coordinates and 
(pt,Pr,P9,P4>) are the corresponding conjugate momenta. 

Next, we can pull the partial derivatives df/dE etc. 
outside of the integral. It is then easy to see that the first 
two terms vanish, since there do exist a conserved energy 
and a conserved z-component of angular momentum for 
the perturbed system. Thus, Eq. (I5.44[) reduces to 



0. 



(5.45) 



Since M(0) is functionally independent of E and L z , the 
prefactor df /dK must be nonzero, so we obtain 



£ e -M(A) = CfjM(X) = 0, 



{K ,Hi} = Q. 



(5.46) 
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The result (|5.46|) applies to fully relativistic orbits in 
Schwarzschild. We need to take the Newtonian limit of 
this result in order to use the result we derived in the 
previous section. However, the Newtonian limit is a lit- 
tle subtle since Newtonian orbits are closed and generic 
relativistic orbits are not closed. We now discuss how the 
limit is taken. 

The integral (|5.46[) is taken over any closed orbit in 
the four dimensional phase space (r,6,p r ,pg) which cor- 
responds to a geodesic in Schwarzschild. Such orbits are 
non generic; they are the orbits for which the ratio be- 
tween the radial and angular frequencies u> r and tog is a 
rational number. We denote by q r and qg the angle vari- 
ables corresponding to the r and 9 motions [36[. These 
variables evolve with proper time r according to 



qr = ?r,0 + kVT, 

Qe = qex> + ueT, 



(5.47a) 
(5.47b) 



where q r ,Q and qg$ are the initial values. We denote the 
integrand in Eq. (|5.46l) by 

l{q r ,q e ,a,e,L), 

where X is some function, and a, e and t are the parame- 
ters of the geodesic defined by Hughes [HJ (functions of 
E, L z and K ). The result (|5.46p can be written as 



T/2 



T/2 



dTl[q r (T),q e (T),a,e,L] = 0, (5.4 



where T = T(a, e, u) is the period of the r, 6 motion. 

Since the variables q r and qg are periodic with period 
27r, we can express the function X as a Fourier series 

oo 

l(q r ,q e ,a,e,L)= £ l nm (a, e, [)e^ . (5.49) 



n.rn— — oo 



Now combining Eqs. pTiTJ) . ([5^5]) and ([5TI9"]) gives 



= J2 Znm(a,e,L)e m ^° 



,o+imqg,o 



n.m— — oo 



xSi [{nuj r + mujg)T/2] 



(5.50) 



where Si(x) = sin(x)/x. Since the initial conditions q r ^ 
and qgfi are arbitrary, it follows that 

X nm (a, e, t)Si [(nui r + mujg)T/2] = (5.51) 

for all n, m. 

Next, for closed orbits the ratio of the frequencies must 
be a rational number, so 



w r p 
wg q ' 



(5.52) 



where p and q are integers with no factor in common. 
These integers depend on a, e and i. The period T is 



given by 2tt/T = qco r — pug. The second factor in Eq. 
(|5.51[) now simplifies to 



Si 



(np + mq)ix 



pq 



(5.53) 



which vanishes if and only if 



n = nq, m = mp, n + m/ 0, (5.54) 
for integers n, fa. It follows that 



Inm(a, e, l) = 



(5.55) 



for all n, m except for values of n, m which satisfy the 
condition (I5.54[) 

Consider now the Newtonian limit, which is the limit 
a — > oo while keeping fixed e and i and the mass of the 
black hole. We denote by I^(q r , qg, a, e, t) the Newtonian 
limit of the function I(q r ,qg, a, e, t). The integral (|5.48[) 
in the Newtonian limit is given by the above computation 
with p = q = 1, since oj r = cog in this limit. This gives 

. oo 

-fdTl N = ZNn,-n(o,e,Oe in( *- - w ' o) , (5.56) 

n— — oo 

where Xn n m are the Fourier components of Xn • In the 
previous subsection, we showed that this function is non- 
zero, which implies that there exists a value k of n for 
which X Nfe! _ fe 7^ 0. 

Now as a —¥ oo, we have uj r /ujg — > 1, and hence from 
Eq. (|5.52[) there exists a critical value a c of a such that 
the values of p and q exceed k for all closed orbits with 
a > a c . (We are keeping fixed the values of e and i). It 
follows from Eqs. (|5.54p and (15.55)) that 



Xfc,_fc(q,g, t) 

XNfc,-fc(a,£,0 







(5.57) 



for all such values of a. However this contradicts the fact 
that 



X 



Nk,-k 



(a,e,i) 



(5.58) 



as a — > oo. This completes the proof. 

Hence, if the three assumptions listed at the start of 
this subsection are satisfied, then the conjecture that all 
vacuum, axisymmetric spacetimes possess a third con- 
stant of the motion is false. 

Finally, it is sometimes claimed in the classical dynam- 
ics literature that perturbation theory is not a sufficiently 
powerful tool to assess whether the integrability of a sys- 
tem is preserved under deformations. An exam ple that 
is often quoted is the Toda lattice Hamiltonian [38l [39| . 
This system is integrable and admits a full set of con- 
stants of motion in involution. However, if one approx- 
imates the Hamiltonian by Taylor expanding the poten- 
tial about the origin to third order, one obtains a sys- 
tem which is not integrable. This would seem to indicate 
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that perturbation theory can indicate a non-integrability, 
while the exact system is still intcgrable. 

In fact, the Toda lattice example does not invalidate 
the method of proof we use here. If we write the Toda 
lattice Hamiltonian as H(q, p), then the situation is that 
-ff (Aq, p) is integrable for A = 1, but it is not integrable 
for < A < 1. Expanding H(Xq, p) to third order in A 
gives a non-integrable Hamiltonian. Thus, the perturba- 
tive result is not in disagreement with the exact result 
for < A < 1, it only disagrees with the exact result for 
A = 1. In other words, the example shows that pertur- 
bation theory can fail to yield the correct result for finite 
values of A, but there is no indication that it fails in ar- 
bitrarily small neighborhoods of A = 0. Our application 
is qualitatively different from the Toda lattice example 
since we have a one parameter family of Hamiltonians 
H(X) which by assumption are integrable for all values 
of A. 



VI. CONCLUSION 

We have examined the effect of an axisymmetric 
quadrupole moment Q of a central body on test parti- 
cle inspirals, to linear order in Q, to the leading post- 
Newtonian order, and to linear order in the mass ratio. 
Our analysis shows that a natural generalization of the 
Carter constant can be defined for the quadrupole inter- 
action. We have also analyzed the leading order spin self- 
interaction effect due to the scattering of the radiation off 
the spacetime curvature due to the spin. Combining the 
effects of the quadrupole and the leading order effects 
linear and quadratic in the spin, we have obtained ex- 
pressions for the instantaneous as well as time-averaged 
evolution of the constants of motion for generic orbits un- 
der gravitational radiation reaction, complete at 0(a 2 e ). 
We have also shown that for a single multipole interaction 
other than Q or spin, in our approximations, a Carter- 
type constant docs not exist. With mild additional as- 
sumptions, this result can be extended to exact space- 
times and falsifies the conjecture that all axisymmetric 
vacuum spacetimes possess a third constant of motion for 
geodesic motion. 
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Appendix A: Time variation of quadrupole: order of 
magnitude estimates 

In this appendix, we give an estimate of the timescale 
T OV oi for the quadrupole to change. The analysis in the 
body of this paper is valid only when T cvo i 3> T rr , where 



T rr is the radiation reaction time, since we have neglected 
the time evolution of the quadrupole. We distinguish be- 
tween two cases: (i) when the central body is exactly non- 
spinning but has a quadrupole, and (ii) when the central 
body has finite spin in addition to the quadrupole. 

1. Estimate of the scaling for the nonspinning case 

For the purpose of a crude estimate, the relevant in- 
teraction is the tidal interaction with energy 

QijSij ~-^S-Qi cos 2 e, (Ai) 

where Sij is the tidal field, 9 is the angle between the 
symmetry axis and the normal to the orbital plane of 
77i2, and we have written the quadrupole as Q ~ QI, 
where Q is dimensionless and I is the moment of inertia. 
For small deviations from equilibrium, the relevant piece 
of the Lagrangian is schematically 

L~Iip 2 + QI^ 2 . (A2) 

We define the evolution timescale T evo i to be the time 
it takes for the angle to change by an amount of order 
unity, and since the amplitude of the oscillation scales 
roughly as ~ m^/mi, the evolution time scales as 

^-sffl (A3) 

where uj 2 rhit = M/r 3 . Thus, the ratio of the evolution 
timescale compared to the radiation reaction timescale 
scales as 




2. Estimate of the scaling for the spinning case 

When the body is spinning the effect of the tidal cou- 
pling is to cause a precession. For the purpose of this 
estimate, we calculate the torque on mi due to the com- 
panion's Newtonian field. The torque N scales as 

Ni i~ £imjQmk£jk ■ (A5) 

We assume that the precession is slow, i.e. 

w P rec < S/mi (jtt) > (A6) 

where w prec is the precession frequency and S = S/m 2 
is the dimensionless spin. This gives the approximate 
scaling of the precession timescale as (cf. 37]) 

WT rr ~ | (j) ■ (A7) 
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and the evolution timescale is thus 



if we define a new time parameter t by 



T IT 

-*- cvol/ -*■ rr >=; I I ■ 

m 2 Q \ r J 



(A8) 



Because of our assumption (|A6[) that the precession is 
slow, equation (| A8|) is valid only when 



!>(£)? 
VM/ Q 



T* \ 3 

~Mj 



(A9) 



When 5 is sufficiently small that the condition (| A9|) is 
violated, the relevant timescale is instead given by Eq. 

E9. 



3. Application to Kerr inspirals 

For Kerr inspirals, 

S ~ a, Q - a 2 , /i/M < f and r - M. (AfO) 



Therefore, the condition (IA9|) is satisfied, and the pre- 
cession time is longer than the radiation reaction time 
by 



T IT ~ i f— 

prcc/ ■*- rr I 



(All) 



Note that for Kerr inspirals, since r ~ M both formulas 
(|A3|) and (|A7p give the same scaling. 

Moreover, for Kerr inspirals, the amplitude of the pre- 
cession will be small, of order the mass ratio fi/M. This is 
because of angular momentum conservation: in the rela- 
tivistic regime, the orbital angular momentum is a factor 
of n/M smaller than the angular momentum of the black 
hole and can therefore not cause a large precession ampli- 
tude. Even if the orbital angular momentum at infinity 
is large, most of it will be radiated away as outgoing 
gravitational waves during the earlier phase of the inspi- 
ral. This factor of fi/M is taken into account when we 
consider the evolution timescale, which for Kerr inspirals 
reduces to 

Wr„~ (£) (i) (f ) ■ <-> 

Since 1/a > 1, M/r ~ 1 and M/fi <C 1, the evolution 
time is long compared to the radiation reaction time and 
we can neglect the time variation of the quadrupole at 
leading order. 



di = ^rdl 



(Bl) 



This is the analog of the Mino time parameter for 
geodesic motion in Kerr The equations of motion 

(j2~26|) - (|2~24l) then become 



o?f x 2 
di 



Mr), 



Vf(f) = 2Er + 2F - K - 4SL z f 



V 3 0) = K 



Li 



sm 2 d 



QE cos 20, 



dip 
di 



QL Z 



sin 



The parameters t and t are related by: 

dt 
di 



Vtf(r) + V t§ (§) 



(B2) 

(B3) 
(B4) 

(B5) 
(B6) 
(B7) 

(B8) 



Q 



V t ~ e {6) = ^-co S 2Q. (B9) 



Vtf(f) 



It follows from Eqs. (|B2[) and (IB4|) that the functions 
f(i) and 0(f) are periodic; and we denote their periods 
by Af= and As. We define the fiducial motion associated 
with the constants of motion E, L z and K to be the 
motion with the initial conditions f(0) = r m in arid 9(0) = 
Omim where f m j n and 9 m - ln are given by the vanishing of 
the right-hand sides of Eqs. (|B2[) and (|B4|) respectively. 
The functions r(t) and 9(t) associated with this fiducial 
motion are given by 



dr 



±^V f {f) 

® d~9 



= t, 



= t. 



(BIO) 
(Bll) 



±^/V § (0) 
From Eq. ((B8j it follows that 



Appendix B: Computation of time averaged fluxes 

1. Averaging method that parallels fully 
relativistic averaging 

We start by noting that the differential equations 
(I2.26|) and (12.271) governing the f and 9 motions decouple 



t(i) = h+ dt' [VtfW)] + V t §[9(t') 



(B12) 



where to = t(0). Next, we define the constant V to be 
the following average value: 



1 



Af jo 



1 f A » 



r / d t'V tf [f(t')} + — / dt'V tg [9(t% (B13) 



As jo 
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Then we can write t(i) as a sum of a linear term and 
terms that are periodic: 



t(t) = t + ri+6t{t), 



(B14) 



where St(i) denotes the oscillatory terms in Eq. (IB12[) . 

To average a function over the time parameter t, it is 
convenient to parameterize f and 9 in terms of angular 
variables as follows. For the average over 9 we introduce 
the parameter \ by 



cos 2 9(i) = z- cos 2 Xi 



(B15) 



where z_ = cos 2 with z_ being the smaller root of Eq. 
(El: 



z± = h 



K + 3QE ±\]{K- QE) + AQEL 2 2 



(B16) 



and where (3 = 2QE. Then from the definition (|Blip 
of 9 together with Eq. (|B4[) and the requirement that x 
increases monotonically with t we obtain 



dx 
di 



= \/P(z + - z- cos 2 x). 



(B17) 



Then we can write the average over t of a function Fg(t) 
which is periodic with period in terms of x a s 



(Fgh = 



1 



[*(*)] 



where 



Ag Jo y/P {z+ - z_ cos 2 x) 



dx 



sjfi{z+- z- cos 2 x) 



, (B18) 



(B19) 



Similarly, to average a function i<f (f) that is periodic with 
period Af, we introduce a parameter £ via 



1 + e cos £ ' 



(B20) 



where the parameter £ varies from to 2-7T as f goes 
through a complete cycle. Then, 



^f[r(0] 



1/2 



pe | sin£ 



(B21) 



(B22) 



(1 + ecos£) 

The average over t of Ff (t) can then be computed from 
/ 27r ^F f /P(0 



Now, a generic function F~ §[?(t),6(t)] will be biperiodic 

mt: F. s [f(i + Af)J(i + A § )}=F.g[f{t),6{i)]. Combin- 
ing the results (|B18|) and (IB23|) we can write its average 
as a double integral over x and £ as 



A a A? 



dx / rfe 



^#(6,0(x)] 



y/P(z+-z- cos 2 x)P(e) 
(B24) 

To compute the time average of E, L z , and K, we need 
to convert the average of a function over t calculated from 
(|B24|) to the average over t. As explained in detail in 
[9j, in the adiabatic limit we can choose a time interval 
At which is long compared to the orbital timescale but 
short compared to the radiation reaction time. From 
Eq. (|B12p we have At = Tt + osc. terms. The oscillatory 
terms will be bounded and will therefore be negligible in 
the adiabatic limit, so we have to a good approximation 



(E)t = ~{EV t ) { , 



(B25) 



where V t = V t f + V t §, cf. Eq. (|F38|) . and similarly for L z 
and K . 

The explicit results we obtain using this method are 
given in section [TTTJ Eqs. 1(535)1 . (|535)l , and (|530)l . 



2. Averaging method using the explicit 
parameterization of Newtonian orbits 

To perform the time-averaging using this method, we 
define a parameter £ via 



P 



1 + e cos £ ' 



(B26) 



where the parameter £ varies from to 2ir as f goes 
through a complete cycle. Note that 9 appears in Eqs. 
p,16[) - (|3 . 1 8[) only in terms that are linear in Q, so we 
can write 9 in terms of £ using the Newtonian relation 



x 3 



r cos ( 



•sintsin(£ + f ) 



(B27) 



(B23) that 



Here, £o is the angle between the direction of the peri- 
helion and the intersection of the orbital and equatorial 
plane. Similarly, for the f9 terms in Eqs. (|3.17p and 
(I3.26[) we can use the Newtonian relations r = e/^/psinf 
and £ = ^/p/r 2 . ^From Eqs. (|2~27|) and ([B20| it follows 
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dt 



o3/2 



d£ (l + ecos£) 2 
and from Eq. (|2~T2|) 



|l - ^2 [-3 + e 2 - 2ecos£ + 2 cos 2 t(8 - e 2 + 8ecos£ + e 2 cos2£)] | , (B28) 



dt 
dt 



Q_ 

2p 2 



(1 



e cos 



[2sin 2 t sin 2 (e + eo)-l]| 



Using these expressions, we compute the timc-avcragccl fluxes from 



and obtain: 



32 (1 - e 2 ) 3 / 2 



(E) 



d£,E(dt/dt){dt/d£) 



J* n dt(dt/dt)(di/dO 



(B29) 



(B30) 



73 



37 4 5/73 823 , 949 4 

1 H e 2 H e 4 ^ 1 e 2 H e 4 

24 96 p 3 / 2 \12 



24 



32 



491 
192* 



P 
P 2 



Q 2 

~2 e 

p 2 



1 85 2 349 4 107 
- H e H e H e 

2 32 128 384 

13 247 2 299 4 



(> 



cos(i) 
179 



11 273 , 847 A 

1 e e 4 H e v 

4 16 64 192 



192 384 
869 1595 



512 



39 
1024' 



e 6 - 



1 



19 



192 384 



23 
512 



e 4 + 



cos(2t) 
3 



1024 



cos(2t) 



96 



48 
1 ; 5 
384 + 384 ( 



121 

— e 4 ) cos(2£ ) sin 2 i 



2084 



e I cos(2£ ) sin i 



(B31) 



32 (1 - e 2 ) 3 / 2 



9 7/2 



COS £ 



l + -e 

g 
16p 2 

5 2 
16p 2 



7 2 5 [61 2 271 4 

8 6 2p 3 / 2 cos t \24 + + 64 6 



61 91 2 461 4 , 
— + — e 2 + — e 4 cos(2i) 



64 



| -3 - ^e 2 + ^e 4 + (^45 + 148e 2 + ^-e^j cos(2t) 

,34] Q , , , 2 /201 51 
l + 3e 2 + - e 4 --^cos^sin 2 , — + -e 



(B32) 



64(1 -e 2 ) 3 / 2 



5 /97 



2p 3 / 2 



37e 2 



1 



55 



139 4 



+ - 



P 

5 2 r 13 



e" H e* + 

2 48 192 



211 
13 841 



cos(t) 



2 , 449 4 

1 e -\ e 4 

4 96 192 



P 



13 



1 



1 



192 64 



512 



V 192 64 



1 

512 e 



cos(2i) 
cos(2i 



Q /391 37 ,\ , fnt . . . o 
-^(l8- + 24 e j 6 C ° S » Sin ^ 



(B33) 



In the adiabatic limit, the terms involving cos(2£o) can 
be omitted because they average to zero. As explained 
by Ryan [l5l |. the radiation reaction timescale for terms 
involving £ is much longer than the precession timescale 
for most orbits, so the terms involving £ W1 H average 
away. This is consistent with our results for the adia- 



batic infinite time-averaged fluxes using the Mino time 
parameter. The Mino-time averaging method was based 
on the assumption that the fundamental frequencies are 
incommensurate and the motion fills up the whole torus, 
which is equivalent to averaging over £ . 
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